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Abstract 

An all electronic 3-D Digital Particle Image 
Velocimetry (DPIV) system has been developed for use 
in high velocity (supersonic) flows. Two high 
resolution CCD cameras mounted in a stereo viewing 
configuration are used to determine the out-of-plane 
velocity component from the difference of the in-plane 
velocity measurements. Double exposure image frames 
are acquired and Fuzzy inference techniques are used to 
maximize the validity of the velocity estimates obtained 
from the auto-correlation analysis. The CCD cameras 
are tilted relative to their respective lens axes to satisfy 
Scheimpflug's condition. Tilting the camera film plane 
ensures that the entire image plane is in focus. 
Perspective distortion still results, but can be corrected 
by proper calibration of the optical system. A 
calibration fixture is used to determine the experimental 
setup parameters and to assess the accuracy to which 
the z-plane displacements can be estimated. The details 
of the calibration fixture and procedure are discussed in 
the text. A pair of pulsed Nd:YAG lasers operating at 
532 nm are used to illuminate the seeded flow from a 
convergent nozzle operated in an underexpanded 
condition. The light sheet was oriented perpendicular 
to the nozzle flow, yielding planar cross-sections of the 
3-component velocity field at several axial stations. 
The key features of the supersonic jet are readily 
observed in the cross-plane vector plots. 


Introduction 

PIV is a widely used technique for measuring the in- 
plane two-component velocity field across a planar 
region of a flow field illuminated by a pulsed laser light 
sheet. Particles entrained in the fluid scatter light 
which is collected normal to the plane of the light sheet. 
The collected light is imaged onto a photographic plate 
or a CCD camera. Multiple images of the particles are 
recorded via short time intervals between laser light 
sheet pulses. The exposures may be recorded either on 
a single photographic plate or on a CCD image frame. 
When CCD cameras are used to record the collected 
light, single exposures may also be recorded on separate 


image frames enabling more flexibility in the 
subsequent data processing. 

Multiple image frame recording can be processed via 
either cross-correlation, particle tracking or a 

combination of both techniques depending on the seed 
particle density. For high seed density flows correlation 
techniques are typically employed where the average 
velocity across a small subregion is computed. 1,2 This 
process is repeated over a regular grid across the image 
frame yielding the velocity vector map. In low seed 
density cases, particle tracking techniques are employed 
to determine the velocity of individual particles on the 
image frame. 3,4 In some instances, both techniques can 
be used to improve on the spatial resolution of the 
resulting velocity vector map. 5 6 

In cases with multiple exposures recorded on a single 
image frame the data processing must be done using an 
auto-correlation technique. Again the average 
velocities across small sub-regions of the flow are 
determined over a regular grid. The auto-correlation 
operation places certain restrictions on the minimum 
displacement of the recorded particle images between 
exposures, therefore reducing the dynamic range. 1 

Although PIV has made 2-component planar 
measurements rather commonplace, many fluid 
engineering problems of interest are fully three 
dimensional and therefore do not possess a plane of 
symmetry for easy 2-component, planar measurements. 
Many approaches have been proposed and investigated 
for obtaining 3-component velocity information 
utilizing the PIV technique. These approaches can be 
broken into two separate classes: 3-component velocity 
measurements in a volume; and 3 -component velocity 
measurements across a plane. Most of the volumetric 
approaches are limited to tracking a relatively low 
number of particles in the volume and suffer from 
particle overlap ambiguities. 7 ' 9 Holographic techniques 
have also been proposed and demonstrated for volume 
3-component velocity estimates via particle tracking. 10 
Barnhart, Adrian ,et. al. have demonstrated a 
holographic PIV system employing high seed particle 
density and correlation data reduction. 11 
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The 3-component planar measurement region 
techniques are simply an extension of the in-plane 2- 
component PIV technique. In the planar approach, 
multiple recording cameras in a stereo viewing 
configuration or creative techniques for combining 
multiple views on a single camera are used. 12 ' 14 There 
are two primary configurations for stereo PIV: lens 
translation; and angular displacement, as shown in 
figure la and b. Lens translation is simple and for small 
field angles results in no appreciable optical distortion, 
but the common image area from the two cameras may 
be restrictively small unless the field angles are 
enlarged. 13,15,16 The use of large field angles in the 
displacement technique increases the sensitivity to the 
out-of-plane displacement, but will also produce 
perspective distortion in the image. 17 

The angular displacement technique is considered to be 
more accurate; however, this technique results in 
perspective distortion on the image plane and typically 
requires a large depth of field to avoid defocusing 
across the image plane. 10,13 The image plane of best 
focus does not fall on the straight back camera film 
plane, resulting in blurring at the edges of the image, 
see figure lb. The use of high f/number collection 
optics increases the depth of field and minimizes the 
blurring. The errors of perspective distortion can be 
ignored by keeping the coupling angle between the 
cameras small. However, sensitivity to the out-of-plane 
component is then reduced. 18,19 A back reflection 
technique has also been proposed to completely remove 
the perspective distortion. 10 

The large depth of field requirement of the angular 
displacement technique is easily eliminated by tilting 
the film plane such that the Scheimpflug condition is 
satisfied. 10,17,20 The Scheimpflug condition requires that 
the image plane, object plane and the median plane 
through the lens all intersect at a common point, see 
figure 2. The tilted object plane is imaged onto a tilted 
image plane where all points are in focus and the depth 
of field requirements of the optical system are reduced. 
The object to image distances vary horizontally across 
the image, resulting in perspective distortion. Both the 
vertical and horizontal magnification vary as a function 
of horizontal position on the image plane. The resulting 
two camera views of a rectangular object are shown in 
figure 2. These distortions can be accounted for in the 
data processing. The advantages of the angular 
displacement technique are the large field of view and 
the improved accuracy in estimating the out-of-plane 
displacement. An indirect advantage of the stereo 
viewing configuration is the collection of the particulate 
scattered light from an angle other than 90° from the 


scattering plane. The scattering efficiency reaches a 
minimum at 90° collection for the case of /^-polarized 
light and small (-1pm) particles, therefore, the oblique 
collection angles collect more light than would be 
detected from a standard PIV setup. 

In this work we will describe and demonstrate an 
angular displacement, stereo viewing, tilt-back camera 
optical system satisfying the Scheimpflug condition. 
The equations relating object to image space are 
derived and used to perform an error analysis on the 
performance of the technique. The calibration 
equipment and procedure will be described. A set of 
validation tests are performed to verify the accuracy of 
the displacement estimates. Auto-correlation PIV data 
reduction is used to process the acquired images and 
Fuzzy inference techniques are used to improve the 
quality of the processed velocity vector maps. 21 
Erroneous velocity vectors are corrected by using the 
Fuzzy inference engine to examine the secondary and 
lower amplitude peaks in the correlation plane and 
replace the erroneous vector with the correct velocity 
vector based on information from the surrounding 
velocity vectors. In this manner the correlation plane 
peak detection is treated as a pseudo particle tracking 
process where all of the peaks in the correlation plane 
are treated as candidate displacement peaks. Flow 
continuity rules govern the Fuzzy inference engine to 
correctly identify the true correlation displacement 
peak. An automatic thresholding operation is also 
demonstrated on the processed data to maximize the 
data quality. Finally, the system will be used to map a 
supersonic flow from an underexpanded jet using cross- 
plane slices through the jet flow. The system is shown 
to perform well even in the presence of the large 
velocity gradient caused by the Mach disk in the 
underexpanded flow. 

Optical System 

An all electronic 3-component PIV system was required 
to satisfy the needs of expedient processing and 
convenient storage of the velocity data. Although high 
resolution CCDs (1320x1040) were to be used, the 
resultant field of view is still quite small when 
compared to that available from photographic plates. 
Therefore, the angular displacement method was 
selected for the optical configuration so as to maximize 
the field of view. The two drawbacks of the angular 
displacement technique are that a large depth of field is 
required to keep the object plane in focus and the image 
magnification varies with horizontal position on the 
image plane. The resulting image of a rectangular 
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object is an isosceles trapezoid, see figure 2. By 
invoking the Scheimpflug condition, a technique 
commonly used in photography for image rectification, 
the large depth of field requirement is eliminated. The 
tilted object plane is exactly imaged onto the tilted 
image plane. Perspective distortion still remains in the 
image. 


The amount of film plane tilt required varies depending 
on the object ( d 0 ) and image (dp distances and the 
angular displacement half angle p. The total coupling 
angle between the cameras is 2 p. The film plane tilt 
angle a is given by: 


a = tan 


d. tan p 

*7 


(i) 


The transformation from object to image space can be 
readily developed by invoking the colinearity condition 
which requires that a ray connecting an object (x 0 ,y 0 ) 
and image point (Xi,yj) passes through the center of the 
lens as shown in figure 2. The object space coordinate 
system is centered at the intersection point of the two 
optical systems at the z o =0 plane. The y-direction is out 
of the plane of the paper. Each camera film plane has a 
coordinate system centered about the middle of the film 
plane at ( x fp ;’y»X Again on the film plane the y- 

coordinate is out of the plane of the page and x lies in 
the plane. The resulting transformations for each 
camera view are: 


d.[x 0 cosPj±z 0 smPj) 


| d o -(Z'Cosfij + jr 4> sin/J > )±tancr x 0 cos/? y ±z a sin/3,]} 


(2) 


d,y 0 


' ,J | d 0f - ( z a cos P j + x 0 sin P y ) ± tan a j [x o cos p } ± z 0 sin 0,]] 




d[(x 0 +dx)cosp l +(z tf +dz)sm^ l ] 

x t +dx t = — +x (5) 

d, ~(z 0 +dz)A t +(x l ,+dx)C t 


where C]=sin pi + tan ( X \ cos P\ y 

A]=cos p\ - tan oq sin Pi and 

dxij is the perceived displacement on the 

camera sensor surface. 


Now by subtracting equation 4 from equation 5 and 
assuming terms of order dz 2 , dx 2 , dzdx « d o we find: 


dx i = M , jctc cos p x 
dzsinp l 


/ 0 v Af,C. 

1 ~ + z o tan P x ) “ 


1 + 


z +- 


tan/3, 


MA 


( 6 ) 


M.A. M.C. 

dy,=MAdy+dzy— dxy o ——\ (7) 

a d. 


and similarly for camera 2: 


dx.=M, 


dx cos P 2 


dz sin p. 


( v MjC. , 

1 + (jc. - z, tan/J, )-*-*• 


1+ « 

l, tan p 2j 


M A 


(8) 


where 7 = 1,2 corresponds to the two camera views, and 
the top symbol in the ± corresponds to view #1 and the 
bottom symbol to view #2. A more useful form of the 
equations would be to express the object plane 
displacements in terms of the two measured image plane 
displacements (one set from each camera). The next 
few steps are designed to achieve the desired result. If 
we assume the particle moves to the new object plane 
location x 0 +dx ; y 0 +dy , Zo+dz then on the image plane 
we can write for camera #1 : 

x T (4) 
d or z o\ +x o c i 


M,A, , MX, 
dy, -MAdy+dz-y o +dx y 0 

l d H d ‘, . 


(9) 


where C 2 =sin fh + lan a 2 cos fo, 

A 2 =cos $2 - tan a 2 sin and the optical 
system magnifications are: 


and 


M,= 


d o l -^o A t +X o C I 


m 2 = 


- Z A - X o C 2 
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These equations describing the recorded particle 
displacements on the image planes can be inverted to 
get the explicit expressions describing the particle 
displacements in the vicinity of the object plane Ax, Ay, 
A z: 


A z = 


sinft 


1 + 


dx x dx 2 
A/, M 2 1 



M,A, 

+ sin 

1 + 

f "S 

x 

z 0 — 


{ t » aJ 

1 d. 

*1 



tan^J 

I d 

4 


( 10 ) 


and d x and d D from both camera views are identical) the 
errors in A z, Ax and Ay are given below: 
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These equations describe the actual particle 
displacements in the vicinity of the object plane in 
terms of the measured displacements on the tilted CCD 
sensor surfaces in cameras 1 and 2 and the location of 
the particle in the object space. Equations 2 and 3 
above can also be inverted so that the two stereo camera 
image coordinates (x, ,y. ) (J=\,2) of a selected object 

plane point (x o ,y 0 ,z o ) can be calculated. 


Error Analysis 
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The equations above describing the physical particle 
displacements in terms of the two in-plane PIV 
displacements can be used to determine the errors in the 
measured velocity components. The computed 
displacements depend on many physical parameters in 
the optical setup which will be ignored in the analysis, 
since their values can be determined to reasonable 
accuracy and their effect should be minimal. Assuming 
that both optical systems are identical (/3 ; = fa, aj -a 2 . 


where 5^;, 5^; are the errors in determining the 

correlation peak locations, and S x , Sy , S z are the errors in 
the location of the correlation subregion on the object 
plane. Some of the partial derivatives have not been 
explicitly shown due to their unruly size. The errors for 
Ax, Ay and A z are plotted in figure 3 over a range of 
coupling half angles for the case y o =0 and Zo=0. The 
independent variable is x Q . The errors in determining 
the correlation peak location were assumed to be 
5dxj=&dyj=Q&5£plMj, where £p is the camera pixel size. 
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Recall that the magnification is not constant in the 
image, hence the error will in general be a function of x 0 
and Zo • Two estimates of Ax and Ay are available, one 
from each camera view. Only one of the two 
expressions for the error in Ax and Ay have been 
included above. In the computations of Ax and Ay the 
average of the two estimates from each camera view 
will be used. 


Figure 3a shows that the error in Ax changes very little 
with horizontal position, growing larger at the 
boundaries of the image and reaching a minimum at the 
center of the image. As the error coupling half-angle 
is increased, the error in Ax increases. Figure 3b shows 
the error in Ay, which is relatively insensitive to 
horizontal position on the object plane and the coupling 
half angle. Increasing the coupling half angle causes an 
increase in the error of Ay only at the boundaries of the 
image. The errors remain relatively flat at the center of 
the image. Figure 3c shows the opposite trends from 
the previous two plots. Increasing the coupling half 
angle /J decreases the out-of-plane displacement error in 
<j aZ . The error in A z is relatively insensitive to 
horizontal position on the image. 

The major variables affecting the out-of-plane 
displacement A z in equation 10 above can be 
determined by making some further simplifying 
assumptions. In addition to the assumptions made 
above we can also assume that 5 ^ = 5^ and x o =0 then: 


a 


Az 


V2Msin /3 


(16) 


where M=d/d a . Hence, we obtain the obvious result 
that the error in the out-of-plane displacement is directly 
related to the error in estimating the correlation peak 
location on the CCD sensor and inversely related to the 
product of the sine of the optical system coupling half 
angle and the optical system magnification. Simplifying 
equation 1 1 in a similar form and combining with 
equation 16 we find: 


a 


Az 



(17) 


This simplification allows us to define an approximate 
ratio of the error in Az to the error in Ax. For the case 
where /?= 30°, the error in Az is 1.7 times as large as the 
in-plane error a AX . At the limit of /3=45°, the error in 
the out-of-plane displacement equals the in-plane 
displacement error. 


In the preceding analysis we assumed that there was no 
limitation on the minimum displacement that could be 
detected by each of the stereo viewing cameras. In 
reality only the cross-correlation technique possesses 
this property. Auto-correlation requires that the 
particles displace by at least their diameters so that the 
displacement peak can be distinguished from the zero- 
order peak. If the flow velocity vector happens to point 
directly along the optical axis of one of the cameras, the 
auto-correlation technique will yield a “dead-band” 
region where particle displacements can not be 
detected. The “dead-band” feature of the auto- 
correlation technique has been ignored in the error 
analysis. The experiments below were designed so that 
the auto-correlation analysis would detect measurable 
displacements. 


Calibration System/Optical Hardware 


The optical system is shown in figure 4. The cameras 
are mounted on carts with sliding dogs that fit into 
concentric radius grooves in the mounting table. The 
center location of the radial grooves is located 8 inches 
from the edge of the table. This configuration permits 
easy changes in the optical system coupling angle while 
maintaining a common focal point between the cameras. 
The lenses are mounted directly on the sliding carts and 
are directed normal to the radial grooves in the table. 
The cameras bodies are mounted independent of the 
lens on the sliding carts and are on rotatable bases so 
that the proper tilt angle a can be introduced. A 
flexible baffle is placed between the lens and the 
camera to prevent stray light from reaching the CCD 
sensor. 

A calibration target consisting of an optical flat with a 
reticle grid on the front surface is mounted in a jig that 
connects to the camera mounting table to hold the target 
at center of rotation of the optical system. The reticle 
grid contains 5 horizontal and 5 vertical lines on 5 mm 
centers. The reticle is mounted on a 2-axis micrometer 
stage on the calibration jig support. The stages allow 
the calibration target to be moved to multiple object z<>- 
planes. 

The optical system is setup so that the middle vertical 
and horizontal lines on the reticle are at the center of the 
camera system field of view. The images from the two 
cameras are electronically subtracted to aid in fine 
tuning the alignment. The center vertical and horizontal 
lines are made to overlap as best as possible with the 
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calibration target at the Zo=0 plane. Calibration data are 
then taken at multiple ^-planes (0, ±0.5, ±1 mm). 

The calibration images are processed and used to 
determine the overall optical system setup parameters. 
The images are first processed using edge detection 
filters to find the horizontal and vertical lines on the 
ruled grid. The slopes and intercepts of all the lines 
(horizontal and vertical) are then determined in a global 
fitting operation. Knowledge about the vertical lines 
being parallel to each other and also perpendicular to 
the center horizontal line are incorporated into the 
fitting procedure. The slope of the center horizontal 
line is used to determine and later correct for the 
amount of rotation in each camera image about the 
optical axis. The extrapolated horizontal lines intersect 
at a point in space which is used to define a keystone 
transformation back to a Cartesian object space using 
just the slopes of either the top or bottom horizontal 
lines. 22 The simple geometric keystone transformation 
is used to determine the optical system magnification at 
the center of the object plane x o =0, y o =0, Zo= 0. 

The intersection points of all of the lines in the 
calibration target for each camera view at all of the z<>- 
planes are then used along with the optical system 
magnification (for each view) to fit for the other optical 
system parameters: a 9 p,d t ,d and the center 

coordinates of the object space on the CCD sensor 
plane . A function minimization procedure is 

used to minimize the distances between the 
reconstructed object plane coordinates of the 
intersecting lines on the calibration target at the multiple 
z-planes and the known 3-D coordinates of the 
intersecting lines on calibration target at the 
mechanically displaced z-plane locations. The fit 
parameters constitute the stereo view optical system 
calibration allowing for the transformation from image 
to object space, or vice versa. 


Equipment and Data Processing 

A pair of Nd:YAG lasers with 200 mJ/pulse energy at 
532nm are used to provide light sheet illumination. 
Two Kodak Megaplus 1.4 cameras are connected to 
EPJX framegrabber boards for data acquisition into a 
single PC. The camera resolution is 1320x1040 and the 
pixels are 6.8^im square with a 100% fill factor. 
Double exposure images are acquired with the 
illumination/recording system. The acquired image 
frames are processed via auto-correlation and a Fuzzy 
logic correlation peak detection algorithm is employed. 


The image frames acquired with the angular 
displacement optical viewing system contain 

perspective distortion where a rectangular object 
becomes an isosceles trapezoid. The perspective 

distortion has two main components: an anamorphic 
component; and a shearing component. Anamorphic 
distortion manifests itself in a foreshortening of the x- 
dimension by approximately the cosine of the object 
plane tilt angle. The shearing distortion is caused by the 
variation in the y-plane magnification with horizontal 
position. The shear is largest at the upper and lower 
boundaries of the image and smallest along the 
horizontal centerline. In order to obtain a regular grid 
of velocity measurements on the object plane, the 
location and spacing of the image plane correlation 
subregions must vary with position. A regular grid of 
points is defined on the object plane and the 
corresponding set of image plane points for each 
camera view is computed and used to determine the 
center location of the correlation subregions. The 
computed velocity vectors are distorted by an amount 
dependent on their position on the image frame. On 
average, the anamorphic distortion for the 30° viewing 
angle is approximately 15%, while the shear is 
approximately 3%. For small subregion sizes, the local 
distortion remains relatively constant across the 
subregion. The distortion is acceptable if it is constant 
across the subregion since the distortion is removed 
when the displacement vectors are mapped back to the 
object plane in the 3-D reconstruction process. Also, 
the 3-D reconstruction is initiated by constructing a 
planar grid in the object space from which the 
correlation subregion center locations are computed for 
each camera view. The average displacements across 
the subregions are computed and used to determine the 
out-of-plane displacements. Hence, the reconstruction 
does not rely on accurately matching the object plane 
coordinates of two image plane points (i.e. particle 
images) to estimate the out-of-plane displacement. 
Rather the out- of-plane displacement is computed from 
the difference of the average in-plane displacement over 
an extended subregion in each camera view, which is 
only sensitive to the order of the size of the subregion. 

During the course of the auto-correlation data reduction 
software development, the choice of the image 
threshold level applied prior to auto-correlation 
processing was found to significantly affect the quality 
of the resulting vector map. The threshold level affects 
the amplitude and width of the zero-correlation peak 
and hence affects the ability of the correlation peak 
detection software to discern the first order particle 
displacement peak from the zero order peak. The 
variation in particulate seed across the image resulted in 
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widely varying background noise levels. A high 
threshold level may produce good results in a high noise 
area, but will yield no correlation results in low noise 
regions of the image. An automatic threshold algorithm 
was found to yield good results over both high and low 
noise regions of the image. The auto-threshold 
technique simply computes the ratio of the number of 
non-zero pixels to the total number of pixels in each 
subregion. The threshold level is continuously 
increased until a ratio of 5% is achieved. The auto- 
correlation operation is then performed on the 
thresholded subregion. 

The velocity vector map obtained using auto-correlation 
processing typically results in a few erroneous vectors 
in the velocity vector map. Many factors influence the 
ability of the data reduction software to accurately 
detect the true secondary peak corresponding to the 
average particle displacement across the subregion. 6 
For each subregion processed, the five largest peaks 
(other than the zero-order peak) on the correlation plane 
are detected and stored. After all of the subregions 
have been processed the Fuzzy inference engine is 
applied. At each subregion a set of velocity vectors are 
detected, which are then examined to determine the 
most probable displacement vector for that subregion. 
The velocity vectors in the surrounding subregions are 
also used in the evaluation of the candidate list at a 
given grid point. In this manner the correlation plane 
peak detection is treated as a pseudo particle tracking 
process where all of the peaks in each subregion’s 
correlation plane are treated as candidate displacement 
peaks. Flow continuity rules govern the Fuzzy 
inference engine to correctly identify the true 
correlation displacement peak. 

Validation Tests 

Tests were performed to validate the capability of the 
stereo view optical system to accurately reconstruct the 
recorded three-dimensional displacements. The goal of 
the tests were twofold: first to determine if the optical 
system calibration was accurate; and second to 
determine if there were any residual effects of the 
perspective distortion in the reconstructed object field. 
Rather than use a computer simulation which neglects 
the effects of optical noise, digitization error, and lens 
distortion, the tests were designed to use as much of the 
actual experimental hardware as was feasible. The 
camera system was configured with a coupling half 
angle j3=30°. A card with a simulated field of randomly 
sized and spaced small white specks on a black 
background was mounted on the face of the calibration 
target reticle. The recorded images contain particles 


that are up to 4x4 pixels in size and there are 
approximately 8 to 10 particles per 32x32 subregion. 
The simulated particle field could then be moved via the 
2-axis translation stage so that the Zo and x Q positions of 
the simulated particle field could be altered between 
image recordings. The single exposure images were 
then electronically added to produce a double exposure 
image for each camera view. 

The simulated particle field was moved through a 
variety of Ax and Az displacements. The simulated 
double exposure images were then processed using the 
auto-correlation data processing software. The two 
velocity vector image maps obtained from the stereo 
views were then processed via the 3-D reconstruction 
software to compute the displacement field. Figure 5 
shows the displacement field results from a particulate 
field mechanical displacement of (0.5, 0.0, 0.3 mm). 
The correlation subregions size was 64x64 with 50% 
overlap. The vector field plotted shows the Ax and Ay 
components of the displacement in the magnitude and 
angle of the vectors. The out-of-plane displacement A z 
is encoded in the color of the displacement vectors. 
Some spurious vectors are observed in the displacement 
map which resulted from low particulate seeding in 
these areas. The averages of the measured 
displacements over the entire field are: Ax:0.493 ± 
0.034 mm, Ay: 0.004 ± 0.011 mm, Az:0.299 ± 0.059 
mm. The absolute value of the Ay displacements were 
taken before computing the mean values to ensure that 
any symmetric distortions in the image did not cancel 
out. The results show that the measured displacements 
accurately reproduce the true mechanical displacement 
of the simulated particulate field. No residual effects 
from the image plane distortion are observed in the final 
reconstructed object plane. The subregion size used for 
this validation test are the largest that would be used to 
reduce the image frames, therefore, perspective 
distortion effects should be negligible. Computing the 
ratio of out-of-plane to in-plane errors yields c^/a^ = 
1.7, which agrees with the ratio predicted from the error 
analysis for an optical system coupling half angle of 
j8=30°. 


Nozzle Data 

The stereo view PIV system was next used to measure 
the 3-component velocity field of a nozzle jet flow. The 
nozzle centerline was directed parallel to the bisector of 
the stereo cameras as shown in figure 6. This optical 
configuration allowed the measurement of cross 
sections of the nozzle flow at various downstream 
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locations. Again the camera half angle was set to 
=30°. The Nd:YAG laser beam was formed into a 2x40 
mm sheet, which was directed from above the nozzle to 
a beam trap below. 

The nozzle consisted of a constant radius contraction 
with an exit diameter of 6.35 mm. The flow was seeded 
in the plenum section using a commercial six-jet 
atomizer filled with Rosco’s smoke juice. The resulting 
mean seed particle size was approximately 1 |im. No 
ambient seed was introduced. The nozzle was operated 
at a plenum pressure of 75 psi, yielding an 
underexpanded supersonic flow. 6 The flow from this 
nozzle has previously been characterized and contains a 
Mach disk approximately 1.5 nozzle diameters 
downstream of the nozzle exit plane. The flow 
upstream of the nozzle is approximately 630 m/s and 
downstream of the Mach disk drops to 150 m/s. The 
flow then repeatedly ramps back up and shocks back 
down. Inter-frame exposure times from 0.8 to 1.2 
microseconds were used to record the velocity field. 
The camera lenses were operated at f/10 and the object 
and image distances were 322 and 47 mm. The imaged 
area was approximately 60x40 mm. 

Cross-plane measurements were obtained at five axial 
stations: 1.1, 1.6, 1.9, 2.6, and 3.0 nozzle diameters 
downstream. The images were processed with an auto- 
correlation subregion size of 32x32 with 50% overlap. 
Six exposure sets at each measurement station were 
acquired and the resulting velocity vector maps 
averaged. Figure 7 shows the resulting 3-component, 
planar cross sections of the measured nozzle flow, 
where the vector magnitudes are color coded. The 
vector plot at z/d=l.l shows the supersonic accelerating 
flow. The cross sections at z/d=1.6 and 1.9 show the 
decelerated flow and the presence of the Mach disk. 
The Mach disk is denoted by the low velocity core flow 
containing the light green/yellow velocity vectors which 
is surrounded by the higher velocity (blue coded) outer 
flow. At z/d=2.6 the flow has ramped back up to 
supersonic and become more turbulent. The final 
station at z/d=3.0 shows again a decelerated flow which 
is also more turbulent than the measurements obtained 
closer to the nozzle exit plane. 

Conclusions 

An all electronic, angular displacement, stereo viewing 
3-component PIY system has been described and 
demonstrated. The image plane focus was optimized by 
employing the Scheimpflug condition via a tilted 
camera back. An error analysis showed the expected 


errors in each of the measured displacements as a 
function of horizontal position on the object plane. 
Residual perspective distortion from the oblique 
viewing direction was demonstrated to be negligible for 
small (32x32) correlation subregions. The out-of-plane 
displacement error was shown to be approximately 1.7 
times larger than the in-plane error for an optical system 
coupling half angle of 30°. Fuzzy inference techniques 
were used to analyze the correlation plane results to 
optimize the validity of the velocity vector maps. The 
stereo view PIY system was demonstrated on an 
underexpanded supersonic nozzle flow, where the light 
sheet was perpendicular to the flow direction. The 
system was shown to be sensitive to the large velocity 
gradients in the nozzle flow. 
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Scheimpflug 

Condition 





Figure 1 : Schematic layouts of the two stereo view 

optical configurations: a) lens translation; b) 
angular displacement. The angular 
displacement cameras are shown with straight 
film backs. The plane of best focus does not 
lie on the film plane. 


Figure 2: Optical configuration satisfying 

Scheimpflug’ s condition. The camera 
backs are tilted by an amount a ; and the 
angular displacement half angle is p. A 
rectangular object is shown to be imaged 
into an isosceles trapezoid. 
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Figure 3: The error in Ax, A? and A z are plotted versus the horizontal position on the object plane. The errors are 

computed for a range of stereo camera coupling half angles (/3=5°, 15°, 30°, 45°). The errors are only 
weak functions of horizontal position. The errors in Ac and A z are strong functions of the optical system 
coupling half angle. The error in determining the correlation peak location on the CCD sensor was 
assumed to be S dx j=8 dy j=0 . 05 £p/M p where Ep is the camera pixel size. 



Figure 4: Orientation of optical hardware on the mounting table. The lens is mounted independent of the camera 

on sliding carts. The carts can be adjusted to change the coupling angle while maintaining a common 
focus point between the two cameras. The calibration target is shown on a two-axis translation stage. 
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Mechanical Displacements: Measurements: 

Ax: 0.5 mm Ax: 0.493 ± 0.034 mm 

Ay: 0.0 mm Ay: 0.004 ± 0.011 mm 

Az: 0.3 mm Az: 0.299 ± 0.059 mm 

Figure 5: Validation test results for a mechanical translation of a simulated particle field. The in-plane 

displacements are depicted in the vector magnitude and angle. The out-of-plane displacements are 
encoded in the color of the vectors as shown by the accompanying key. The actual mechanical 
displacements and the average and standard deviation of the measured results across the plane are 
denoted in the figure. No residual effects from the perspective distortion are observed. 



Figure 6: Experimental setup for the nozzle flow measurements. The light sheet is directed from above to a 

beam dump below the nozzle. An exhaust duct is mounted between the cameras to collect the seeded 
air stream. 
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The flow becomes more turbulent as the axial position increases. 
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